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^«0 . Abstract 

cn 

lO ' We give a complete analytical computation of three and two-point loop integrals occurring 

^^ , in heavy-particle theories, involving a velocity change, for arbitrary real values of the external 

,-_ ' masses and residual momenta. 

^ ' 1 Introduction 

Qh! The study of the dynamics and spectroscopy of hadrons containing a heavy quark has been greatly 






simplified and systematized with the introduction of heavy quark effective theory (HQET) [y. 
Heavy-particle theories along similar lines have also been succesfully applied in other, related 
contexts. Thus, in those cases where a chiral approach to the strong interactions of heavy hadrons 
K> ' with light mesons is applicable, a combination of chiral and heavy-quark symmetries leads to 

j^ ' heavy hadron chiral perturbation theory (HHChPT) Q. A heavy-particle expansion has also been 

developed in the chiral-perturbative framework for nucleon-meson interactions, which constitutes 
the so-called heavy baryon chiral perturbation theory (HBChPT) Q. 

In the heavy-quark limit the interaction of a heavy quark, or hadron, with the light degrees of 
freedom cannot change its four- velocity v^. In consequence, v^ becomes a good quantum number 
and, therefore, heavy-particle effective theories of the strong interactions are expressed in terms 
of velocity-dependent fields. Weak interactions, or other external sources, however, can change 
the velocity and/or flavor of a heavy quark or hadron. Strong-interaction corrections to velocity- 
changing interaction vertices then involve loop integrals with two different velocities. 

In this paper we report on a complete analytic computation of three-point loop integrals in- 
volving a velocity change, and two-point loop integrals. We consider a class of one- loop integrals 
occurring in heavy-particle theories, with arbitrary real values for the external masses and residual 
momenta. Since our aim here is mainly methodological, we will not discuss specific phenomeno- 
logical applications. For definiteness, however, we adopt the language of HHChPT in the sequel. 

In the next section, we define the integrals to be studied, establish our notations and conven- 
tions, and discuss the method we use, which involves a combination of the HQET technique and of 
standard methods for computing loop integrals |Q, |5|, p. In section 3 we give technical details about 
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the computation of the scalar three-point integral, state our results and discuss several important 
limits and particular cases and cross-checks. In section 4 we briefly consider the two-point integral, 
which has already been given in the previous literature. In section 5 the vector and second-rank 
tensor integrals are given in terms of form factors. Finally, in section 6 we give some final remarks. 

2 Method. Notation and Conventions 

The loop integrals we consider are of the form, 

X"l-«n = !^f (fa 1 1 fl) 

{2t:Y J ^{2vi-{q + ki)-SMi+ie){2v2-{q + k2)-6M2 + ie){q^-m^ + ie) ^ ^ 

■ A—d p a\ . . . „o„ 



[I-kY ] ^{2v{q + k)-6M + ie){q^-m^ + ie)' 

Here v^, i = 1,2, are the velocities of the external heavy legs, k^ their residual momenta, and 5Mi 
their mass splittings relative to the common heavy mass of the corresponding heavy-flavor/spin 
multiplet. m is the mass of the light particle, which in HHChPT corresponds to a light pseudoscalar 
meson. These integrals are defined in d = 4 — e dimensions, // being the mass scale of dimensional 
regularization. Their degrees of divergence are n + d — 4 for Xa"^ ■"" and n -\- d — 3 for T2"^"'"". 
The factor of 2 in front of v^ corresponds to our normalization of the heavy-particle propagators. 
In this section we will restrict ourselves to the scalar case n = 0. The cases n = 1,2 will be 
considered in detail in section 5. Together with X2,3 we consider also the auxiliary integrals, 

~ iu^~'^ /■ 1 

^^ " "(2^ J '^ '^{{q+piY - Ml + ie) {{q+P2f - M^ + ie) {q^ - m^ + is) ^^^ 

~ iu^~'^ f 1 

"^^ " l2^ i "^ ^ ((g + pY - M2 + ie) (g2 -m'^ + ie) ' (^) 

X3 is convergent in four dimensions, with degree of divergence d — 6. I2 has degree of divergence 
d — 4, diverging logarithmically at d = 4. The relations among external momenta and masses in 
I3 and I3 are, (z = 1, 2) 

pf = Mv^ + fcf , pfK^ > ; Mi = M + ^6M, , M, > 0, (5) 

and similarly for T2 and T2. We remark at this point that T2,3 need not be related to Feynman 
diagrams in any existing physical theory. The similarity of the limit M ^ 00 studied below with 
the heavy-quark limit is purely formal. X2,3 are just intermediate steps in our calculation of T2,3, 
as we discuss next. 

In the limit M — > 00 we have, {k = 1,2) 

1 1 „M 



{q + PkY - Ml + ie M 2vk -{q + kk)- SMf, + ie \ M^ 

Since 5X2^3 /Qttt-^ are convergent for d = 4, equation ^ leads to, 

m .,2d% , ,nf^\ dl2 _ ,, d±2 .. / 1 



+ 77^ • (6) 



M'lz5i + 0[-], 7d-2=M^,+0- . (7) 



dm? dm? \M I ^ dm"^ dm"^ \M 



Therefore, at d = 4 we must have, 

T3 = XaU^o + m2 (j3 - T3|_o) + O (^^) (8) 

T2 = X2U..0 + M (T2 - X2U.0) + O (^) • (9) 



Moreover, using equations @ and defining Aj = (5Mj — 2vj-kj, we can write T3, 23 in terms of Aj. 
Differentiating we obtain, to leading order in 1/M, 



or, equivalently. 



^r""S-^°i^^' '-'''' ""^ 



T3 = M^I^ + Ci(Ai) + O (^^) = m2J3 + C2(A2) + O (^^) . (11) 



Here, the dependence of C1.2 on d,M,fi,m and fi-W2 is understood, but it is shown expHcitly 
that Ci can depend on Ai but not on A2, and the oposite is true for C2. Subtracting the two 
equations (|ll|) term by term, we conclude that Ci = C2 = C{d,M, 12,171, vi-V2) do not depend on 
Ai_2- (Furthermore, (|^) together with (p^) imply that C does not depend on m either.) Thus, 

l3=l3\A,^o=A,+M^(ls-i3U,^o=A,)+o(^-^^ . (12) 

We notice that X3 is straightforward to compute for Aj = by using the HQET method for 
combining denominators (see, e.g., 0). On the other hand, X3 is needed in ([T^ ) only at d = 4, 
and to leading order in M, including logarithmic corrections. Eq. ( p^ will then be the starting 
point for our computation of 23. Equations for I2 analogous to (|ll]) and ( [T^ ) can also be obtained, 
involving two derivatives. We will find it more convenient to use eq. (H) in order to compute T2. 

Scalar integrals can depend on v^, V2 only through tj = wi •f2- If we denote by O the magnitude 
of the three- velocity associated to Vi or Vg in the rest frame of u f + V2 then S, §11.5], 



ILO 



1 + n' « - v^r I- ^ .,0. 



Together with Q, the roots of (vf — avfi) = 0, given by 



1±0 



a± = u; ± Vc^2_i = __ , (14) 

will appear frequently below. For physical values of fj'2; such that (i"! 2)^ = 1, we have uj > 1, 
0<r2<l,0<a_<l<a+. We will always assume these inequalities to hold in what follows. 

Logarithms have a cut along the negative real axis. The log of a product can be split as log(a6) = 
log(a)-|-log(5) if Im(o) and Im(6) have opposite sign, or if a > 0. Similarly, log(a6) = log(a) — log(6) 
if Im(a) and Im(6) have the same sign, or if a > Q. Given a complex number z, we use that 
determination of the argument such that — vr < arg(2;) < vr. In particular, log(l/2;) = — log(2:). We 
use the same definition and conventions as P] for the dilogarithm, which we denote by Li2. 



3 The Scalar Three-Point Integral 

We will now consider in detail the calculation of X3. Our first step is to compute T3 to leading 
order in M. As mentioned above, we only need to evaluate T3 at d = 4. We introduce a standard 
Feynman parametrization of the integrand in (pi). Integrating over d'^q and over the Feynman 
parameter associated with the third propagator in (P), we obtain, 

1 /•! ri-x I 

(47r)2 io io {ypi + XP2Y + y{Ml — p1 — m?) + xiM^ — P2 — "t-^) + m^ — ie ' 

In the limit M — > 00 the polynomial in the denominator is a sum of terms of the form M^x'^y^~'^, 
with p = 0, 1, 2 and < g < p, the term with p = being m? — ie. These are properties we want to 
maintain, in order to be able to take the limit M — > 00 later, retaining only the leading terms in 
M in each coefficient. Thus, we will not make a change of variable y — > 1 — y at this stage, as the 
limits of integration suggest. Following |^], we shift variables according to y ^ y — ax. This shift 
is homogeneous in x, y, so it does not change the order of each term as M — > 00. The parameter 
a is taken to be one of the roots a± of (p2 — api)^ = 0. For M large, we have a+ > 1 > a_ > 0. 
Choosing a = a+ and exchanging the order of integration, we are led to. 




D = (2yipi ■p2 - a+pj) + (Ml - pj - m^) - a^Mf - pj - m^] 
+Piy^ + y{Mi — p\ — m^) +m^ — ie . 

Using equations (|5[) and retaining only leading powers of M in each coefficient, the previous expres- 
sion simplifies considerably. In the limit M ^ 00, a± are given by (0). Performing the integration 
over X we obtain, 

(17) 



with 



X3 = 


' '-^" ' t\y 


1 , /i^\ 1 i-r^^ 1 r+ , 1 , (^ 


647r2 n M2 io ^ 


g'°^U; 64»2 fi A/2 A ''^o'°su 


L = 


n 4f] y-l 


1 - SZ^ a+ 


H = 


2 Ai m^ . 


^ M ' 




1 + 

yo = -^(^A + 5) 


A = i(Ai + A2); 5 = ^(Ai-A2). 



(18) 



In these last two equations we have introduced several notations that will be needed later. In the 
second integral in (|17|) the variable y is 0(1) over the entire domain of integration. Therefore, 
My = 0{M) and the integral is given, to leading order in M, by 

' dy-log{a_ ^ . (19) 



647r2 n M2 7i " y V -y + a+ + I - ie 
This expression does not depend on Ai^2- Therefore, it will cancel when we subtract IsIa-^^o^a 



from T3, and will not contribute to I3 as given in (12). We shall then drop this term from T3 from 



now on. The remaining integral can be re-written as, 



1 1-^2 1 



dy ^ {log{H) - log{K)} + 



(20) 



•" 647r2 n M2 Jo "" G 

We notice that there is no singularity at the zero of the denominator, since the numerator vanishes 
there. We will denote yi±/M and y2±/M the roots of H and K, respectively. They are given by, 



yi± 



-As ± ^Ai - 4m2 + i^ 



From their definition, (0), it is clear that H and K are equal at the zero of G. Defining, 

Zka = Vka - Vo , k = l,2, a = ± , 

the equality of H and K at y = yo/M can be expressed as, 

an identity that will be important below. 

After factorizing H and K and splitting the logs in (pO|), we find, 

02 1 rl 



(21) 

(22) 
(23) 






M 



log 



a_ 



2/ 



J/2+ 

M 



log 



a_ y 



^2- 

M 



(24) 



In order to be able to distribute the integral inside the braces without introducing spurious singu- 
larities, we use (23) to add and subtract the value of each log at the pole. In this way we obtain, 



^ 647r2 n ^P^^ ' 



dy 



y - Vo/M 



log y 



M 



log 



yo_ 

M 



Vka 

M 



(25) 



where the sum runs over k = 1,2 and cr = it. These integrals are already in standard form. 
Evaluating them to leading order in M, we arrive at. 



1 1-J]2 1 
647r2 n AP 



E(-i 



ifc+i 



^'°<7 



vr 



^log — 
2 V A* 



+ log( -^ ) log 



Zka 



log 



^^fe 



logf jlogf 


Zka\ 


u ) \Zka J 


\Zka 



■ (26) 



In order to simplify this result we have explicitly used the relation \og{fi/zko) = — log(^;fco-//u) as 
explained in section |2|. Notice that we have introduced a mass scale /x that makes the arguments of 
the logs dimensionless. The first two terms in (pq ) are independent of k and therefore they cancel 
out in the sum. The third term can be simplified by making use of ([2^), which results in a term of 
the form (1 - f72)y(647r2M2[7) log(a2) log(M///). 

The integral T3 at Aj = is straightforward to compute. Its logarithmic dependence on M 
cancels that of T3, so that we obtain. 



2^3 — 2r3|A,=o 



1 



1-^2 



327r2M2 2Q 



log2(a^)+Iog(a2)log - +^(-1 



^ 



ilog2f£^ 
2 ^ \ ^l 



M-Y^--) +M--y4-) +M-)]] ■ (27) 

V ^ / \ fJ' J V Zka J \ZkaJ \ZkaJ\} 



This is essentially the final result, except for the dimensional regularization pole term, which is 
supplied by Tal^.^o (see eq. (|T2D). Setting n = 0, 6Mj = = kj, j = 1,2, in (|^), we obtain, after 
using the HQET method for combining denominators and integrating over d'^q, 



2^3 1 A,- =0 



/^ 



4-d 



(4^)d/2 



r 3 



1 /"OO \ 

dx / dX 

Jo [A2 (2;2 + (1 _ x)2 + 22;(1 - x)w) + m2 



le 



i3-d/2 



■ (28) 



The innermost integral can be evaluated by changing variable to u = A^. We get. 



1 r (2 - ^] [^ 

2(47r)''/2 V 2)\m 



4~d pi 



dx- 



1 



1 



2(47r)<^/2 V 2)\m) 4J7 



x^ + (1 — xY + 2x(l — x)iv 



log(a+) 



11-^' . ^ /2 . . /^' 

log(a+) - - 7e + log(47r) + log ^ ^ 



647r2 



m^ 



+ Oie) 



(29) 



where 7e is Euler's gamma. 

Thus, finally, using (|l|), {^) and (|9|), we obtain. 



T, 



_ ^ ^ ^ . ^ Inrrri' ^ 1 Inr^fn- 1 1 Vf I'l'^^ 


"1, 2/^^fe<T 


647r2 n ] € ^°S^''+^ ' ^°S la+J 1 2^1 ^) 
-log ^ log — — +log log + 

\ IJ' J V At / V Zk„J \ZkaJ 



(30) 



where we used (23) to rewrite log(//) in (27) as log(7l) = log(//-v/47r) — 7e- Equation (30) is our 
main result. It gives the analytical expression for T3 for general real values of the external masses 
and residual momenta. Below we consider some particular values of the parameters which are 
important in practice, and in which the expression for T3 takes on a simplified form. They can also 



serve as cross-checks of (30). 

We note that T3, given by (ffl) with n = 0, is symmetric under exchange of heavy particles 
1 <-> 2. The general result (30) is not manifestly invariant under Ai ^^ A2, but its symmetry has 
been thoroughly checked numerically. 



3.1 The case ya < and the zero-recoil limit 

An important particular case to consider is the case Vi = V2 or, equivalently, fi = (w = 1). 
X3 at ri = can be computed more directly by differentiating T2 (see section ^. Therefore, its 
calculation from ( |30|) constitutes a cross-check. 

The limit fi ^ 0, however, is difficult to take in (|30|). As indicated by the factor 1/il in that 
equation, the limit results from complicated cancellations among different terms in the sum in (|30|). 
Individually, some of those terms may be large as O — > 0. This situation arises mainly from the 
addition and subtraction of terms of the form log(yo/M — yka/M) in (^), which is necessary when 
2/0 > 0, since in that case there is a singularity in the integration domain. If, however, we restrict 
ourselves to the case yo < 0, then the general expression (pH) takes a simpler form, which makes 
the zero-recoil limit transparent. 

As mentioned at the end of the previous section, T3 is symmetric under exchange of Ai and 
A2. Hence, without loss of generality we can assume 6 > (see eq. (|lq). The case 6 = will be 



considered afterwards as a limiting case). From the definition of yo in (1^) we see that, if 5 > 0, 
then for sufficiently small 0, we will have yo < 0. This is therefore the relevant parameter region 
to consider when fi ^ 0. 

Assuming, then, yo < 0, we can go back to (|2^) and distribute the integral inside the braces 
without adding extra terms. The calculation goes through unchanged, yielding. 



c5>0 
yo<o 



1 1-^^ 

647r2 n 

+ E(-i)' 



log(a+) 



^ + log(a,) + logf^ 



Vka 



Vq 



log- 1_— +l0g 1_— lO; 



Vko 



Vo 



(2VLy, 

; 

V Vka 



\Zka 



■ (31) 



The zero-recoil limit of this expression can be easily obtained. Using the assumption (5 > and 
the relation (valid for i7 = 0) 

Y.{-lfyk. = 25 , 



we find, 



Z-^ 









Vko 



(32) 



It is understood that in this equation we must set il = in the expression (21) for yj^o-- In section 
^ we will obtain this result from X2- 

We give, finally, the result for T3 at zero-recoil when 5 = (Ai = A2 = A). We just take the 
limit of (^2|) as 5 —> to obtain. 



1 [2 ^ hf\ A 



-^— B^, 



(33) 



where ^/'■ = ^/ IS? — Am? + ie. Notice that, in general, we cannot express the difference of logs in 
as the log of the ratio of the arguments , because their imaginary parts have opposite sign. 



3.2 The case y^, = 

When 2/0 = the general expression ( |30|) for T3 is singular. The singularity is avoidable, though, 
so that we can take the limit yo — ^ in (pO) to get. 



!/0=0 



1 1-n^ 

647r2 n 
^647r2 n 



log(a+ 

E(-i) 



+ log(a4 



-log(^— j-log(^— )lo, 



Vka 



(34) 



This expression is valid, in particular, when Ai = = A2, which is the point in parameter space 



we used to "match" T3 and T3. As is easily seen, in that case we recover (29). 



3.3 The case m = 

The case m = is relevant to theories involving massless particles, such as gluons in HQET and 
Goldstone bosons in chiral theories in the limit of massless quarks . The value of X3 at tti = 



can be obtained from the general expression (pO|), and also by direct computation from (|l|) using 
the HQET method for combining denominators. We will consider both approaches in this section. 
Together with the zero-recoil case studied in sections q^ and ^, this is one of our main cross-checks. 
We will now study the limit m ^ of X3 as given in (|30| ) . We notice that when m = one of 
the two possible values cr = it leads to yka = (see eq. (pil) ). In fact, m = and A^ > (resp. 
Afc < 0) implies yk+ = (resp. yk- = 0). Calling ak = -sgn(Afc) = -a'^, so that yk^k / = y^^/^, 
we have z^o-' = ~?/o + ^EO'fc and, therefore. 



loe 



Vo 



log 



0; 



LiJ^ 



^ka'i^ ) \ ^ku'^ ) \ ^ka'i^ 



0, 



so these terms drop from the sum in (|30|). On the other hand, since yo is by definition (|l^) a real 
number, we have the equalities, 

2 2 

Y^(-l)^^log^-y, + iea',) = m0(yo) log(yo) E(-1)M 
fc=i fc=i 

2 2 

E(-l)'''log(-yo + «ecT^)log(yo-^e4) = i'^sgn{yo)log{\yo\)^{-l)''(Tk , 

fc=l A;=l 

6(x) being a step function. With these relations taken into account, from 
1 1 - (72 f 2 ^ 



we find for T3, 



647r2 [7 



-log! ^ llog 



log(a+)+log2(a+) + E(-l 



£{yo) 



|yoi 



+ log 



fc=l 

2/0 \ 



1 



2 \ fi 



ivr^(-yo)log(^ 5:(-l)'=4 

/^ ^ A:=l 



-2i7rlog 



log^l+Lisf^'"'^ 



-2/0 



A* 



0(Ai)0(-A2) . 



+ SiVo) 



,(35) 



(36) 



Notice that the sum in (^) runs over k = 1,2 but not over ak, and that — yo > in ( [3q ) due to 
the step functions. As before, y^o-fc and Zfco-fe in (HI) are given by (|2l|) and (|2^) with m set to zero. 
When 6 = (Ai = A2 = A), equation ( |35| ) becomes 



ei^^T^ {7 logK) - ^ logK) iog( ^ 



02 



2 log (a+) log 



ze 



/^ 



+Li2 



1-n 



Li, 



1 + S7 



(37) 



a result we will explicitly cross-check below. 

We now turn to the calculation of X3 at m? = directly from its definition (|l|) with the 
HQET method for combining denominators. For the sake of brevity, we will skip the details of the 
derivation and quote the final result, which can be written as, 



I3 



1 1-n^ (2, . , 1, , ,, f^-^^\ T- f^-^ 



Li2 



1 1 



n' 



647r2 n 

A + 6 -ie 



,„„ JL±3L) U^l^ii) , U,d^^3L , ,„, 



-Li2 



A-5/n- ie 



nA-6-ie, 

f A-6-ie 
^^'\A-6/n-ze 






'\nA-6-ie 



l + Q 
2 

A-6 



l£ 



fJ- 



-6) 



(38) 



Here, the terms within braces in the second and third hne are to be repeated with 5 replaced by 
—6 as indicated. 
When 6 



0, (pSf ) reduces to (|37|). Thus, in the case m = = 5 we have an analytic cross-check 
of our results. In the more general case 6 ^ 0, we have numerically compared general expression 
(|30|) for small values m, with equations (|35| ) and (|38|). The three expressions for I^I^^q were found 
to agree over a wide range of real values for Ai, A2, and < < 1, again providing a cross-check 
for @). 



4 The Scalar Two-Point Integral 

The scalar two-point integral X2 is given by (^) with n = 0. T2 is a function of m and A = 6M—2v-k. 
The starting point for the calculation of I2 is (^). The computations of both I2 and 'l2\m=o are 
standard. Defining, 

x± = - (-A ± VA2 - 4m2 + ie) , (39) 



from (P) we obtain, 



32tt' 



- + 2 + 



levr^ 



x+ log 



-^ -hx_log — ^ 



(40) 



This result can also be obtained by using the HQET method for combining denominators, which 
yields an equivalent expression in terms of hyper geometric functions. 

In order to compare our result ( |40| ) for X2 with those in the previous literature, we rewrite it 
in terms of x = A/(2?n,), 



T2(A,m) 



A 
32^ 



U.,{tV. 



m^ 



m 
167r2 



•^(x) 



(41) 



with ^(x) = \/x2 — 1 + ie logfx — \/x-' 



1 + ie) — logf x + \/x2 — 1 + iej . The coefficient of the 
dimensional regularization pole vanishes when A = 0. This is due to the fact that the real part of 
the integrand in (|2|) is parity-odd when A = 0. 

Equation (41) agrees with [^ for all real values of A, once their different normalization and 
conventions are taken into account. It agrees with the results from [10, |ll| (see also the second of 
P]) only in the region x > 0, our result being different from theirs over the entire negative semiaxis. 

We consider now X3 at zero recoil. As shown in [|l^], and seen from its definition (|^), T3 at 
il = can be obtained from T2 as. 



(Ai,A2) = -^ — (X2(Ai,m)-J2(A2,m)), and X3 (A, A) 

Ai — A2 f!=0 



_d_ 
dA 



X2(A,m). 



(42) 

Substituting the value of X2(A,m) given by ( |4l| ) in (|4^), we recover our previous results (^2|) and 
), as can be easily checked. 



5 Vector and Tensor Integrals 

In this section we give general expressions for vector and second-rank tensor integrals in terms of 
form factors. We also compare our results to those in the literature, when available. The form 



factors will be expressed in terms of scalar integrals. Of those, X^ and 23 have been given in 
previous sections. We will also need, 



til L]'^ 1 ^ _rr^ (2 (^ 

'" {iTrYJ '^ {q^ - m? + ie) ~ IQir"^ \e ^ "^^m? ^ ^ 

with ?TT, > 0, d = 4 — e, log(7I^) = log(/i^47r) — 7^. Two other scalar integrals appear in the 
evaluation of tensor ones, 

X/ = !^ [d\- L ^ and I2' = ^^ ld%- -t^ (44) 

{2-kY J ^{2vq-A + ie) (2ttY J ^ {2vi-q - Ai+ ie) {2v2-q - ^2 + ie) ^ ' 

Both Ti', and Z2' vanish, as we will now show. The easiest way to see that Zi = is by applying 
the axioms of dimensional regularization [|l2| , §4.1]. We consider Ti{v^,d) as a function of d and 
v^^, momentarily allowing v^v^ > 0, not necessarily equal to 1. Then, we can always shift the 
integration variable so that, 

4-d 



^■'(""■"'"I^tM 



(27r'^) J ^{2vq + ie) 

Let s > 0, and consider Zi'{sv^,d). By factoring s out of the integral we get, Ii'{sv^,d) = 
l/sTi{v^,d), whereas by rescaling the integration variable q^ we find, Ii'{sv^,d) = 1 / s'^Ti [v^ , d) . 
Therefore, we must have l/sTi{v^,d) = 1/s Zi'{v^,d) for all s > and all complex d, excluding 
positive integer values. Since Ii must be an analytic function of d, we conclude that it vanishes 
for v^Vfj, > 0. 

We now turn to X2'. Introducing a Feynman parameter, we can write it as, 

^2' = ^ f'dx fd%, ^ 



{2ttY Jq J ^2V{x)-q- /\{x)+i£' 

with V^{x) = xvi + (1 — x)v2 and A(x) = xAi + (1 — x)A2. For Vi = 1 = ^2 and < x < 1 we 
have y^^(x) > 0. Thus, the inner integral is Ii'{V'^, d) and since Ii' = 0, I2' also vanishes. 

5.1 Vector and tensor two-point integrals 

The vector two-point integral is given by, 

^r(^.',A, ,n) = |t! |<^,____|____ . (45) 

Lorentz invariance dictates that I2 is given in terms of only one form factor, which can be imme- 
diately evaluated by algebraic reduction [^, p, 

T2''(^^", A, m) = F(A, m)v^ , with F(A, m) = v^l2^{v'', A, m) = ^Xi(m) + ^M^. H , (46) 

where the scalar integrals Xi and X2 have been given in (^) and in section ^, respectively. 
The tensor two-point integral is defined as, 

l2'"'{v'^,/\,m) = ^I-—r ld\- "1^^ . -. (47) 
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We will introduce two sets of form factors. First, we define, 

2'2^'^(w°,A,m) = IoiA,m)g'"' + h{A,m)v^'v''. (48) 

Second, we introduce F form factors which can be easily computed in terms of scalar integrals, 

(49) 
(50) 



Fo(A, m) = g^a2^^ = dio + h = m^^A, m) 
Fi{A,m) = v^v^l2^'' = Io + h = - (Ti(m) + Al2{A, 



m 



In fact, Fq(A, m) = Xi' + 77^^X2 (A, ?7i), so here we have used Xi' = 0. Inverting the relation among 
F's and I's we obtain, to lowest order in e = 4 — d. 



/o(A,m) 



K-l 



A [ A^ \ 

-Xi(m)+ — -m^ X2(A,m) 



A / e \ r 1 m ( e 

h{A,m) = -h+_j [Xi(m) + AX2(A,m)J - _ M + - ) X2(A, m) 



(51) 
(52) 



Finally, we substitute the known values of Xi(?7i) and X2(A,?ti). Using the same notation as in 

i), 



/o(A,m) 
A(A,m) 



m3 f/2 ^ ^f]!^' 
3 • IGvr^ I \ e I ^Ti^ 



X X 



2/ 2 ^ 



l)-^(x) 



m 



3 • levr^ 



^ + log ( ^ ) + f ) x(4x2 - 3) + ^x + (4x2 _ i)_^(^) 



(53) 
(54) 



This is the general form for X2 . We do not find agreement with [0| . 

There are two particular cases of interest, in which X2'^'^ can be easily computed directly by using 
the HQET method for combining denominators, thus providing cross-checks for our results. In the 
first place, we consider the case m = 0, A > (the case A < is analogous). A straightforward 
computation using the HQET method yields, 



If 



A3 



m=0 
A>0 



3 • 1287r^ 



hM^ 



+ 



8 



h 



A3 



'"=" o on 2 I 

A>o 3 • 3z7r^ \ e 



+ log' "^ 



A2 



+ 



13 



(55) 



which agree with (^3|) and ( |54D evaluated at m = 0. 
Second, in the case A = we find, 



X2^'^(i;",A = 0,m) 



m 



3-167r 



(^M^ - v'^v'') 



(56) 



again in agreement with the corresponding limit of ( |53| ) and (54). As remarked above, in connection 
with the scalar integral, there is no dimensional regularization pole in this case. 

Furthermore, if we assume < x < 1 {m > A/2 > 0) and expand in powers of x, we recover 
the result given in Q. 
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5.2 Vector three-point integral 

We now turn to the tensor three-point integrals, starting with the vector one, 



J3'^K,r;^,Ai,A2,m) 



Z/i 



(27r)° 



d% 



{2vi-q — Ai + ie) {2v2-q — A2 + ie) {q^ — m? + ie) ' 



(57) 



On the left-hand side we omitted /i and d from the argument for brevity. We define, as before, two 
sets of form factors, 

Jg/^ = hv^ + /2V2 and Fi^2 = ^^1,2 • 2^3 , 



with /j(J7, Ai, A2,'m) and Fj(r2, Ai, A2, w,) related by, 



1^ 



-{l-a^)Fi + {l + n'^)F2 






(1 + n^)Fi - (1 - a^)F2 



The form factors Fj can be expressed in terms of scalar integrals as, 

Fi,2 = -T2(A2,i,m) + ^J3(J^,Ai,A2,m) 

These equations give an explicit expression for T^^. 
At zero recoil we have v^ = V2 = v^^ and. 



T^^ 



Vu ■ I3 



Using (^ we obtain 



n=o 



v'' , with Vy-T" 



'Jv ■ ^3 



(58) 



(59) 



(60) 



^^^ = ^T2(A2,m) + ^T3(f^ = 0,Ai,A2,m) . (61) 



v^ ■ T3" 



Q=o 2 Ai - A2 



[AiT2(Ai,m)-A2T2(A2,m)] 



if Ai / A2 and 



if Ai = A2 = A. This completes our treatment of the vector integral. 



.^o = ^^^(^'"^) + yM^^(^'"^) 



5.3 The tensor three-point integral 

We consider, finally, the tensor integral. 



g^g^ 



^ ^i'2' i' 2, ; (2Tr)d J ^{2vi-q-Ai + ie){2v2-q-A2 + ie){q^-m^ + ie) 



In this case we have two sets of four form-factors each, 



13^" = hlV^iV\ + l22<^^2' + h2vTvl' + /o5' 



{m^,*^} 



^tiu 



with V1V2 = V1V2 + vt^Vi, and. 



{^,,1'}-- 



Fii = v^v^l^^y ; F22 = v!^V2l3f,^ ; Fu = vfv2^l3^^ ; Fq = g^^'lz^y . 



(62) 



(63) 



(64) 
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The F's can be given in terms of the /'s using (pSl). Inverting those relations we obtain, 

d-i{i-ny 



'11 



d-2 I6J74 
d-2 4rj2 



;i - n^fFn + {n^ + 2^0^ + 1 ) F22 - (1 - n')F, 



d-l 



12 



+ 



Fn 



j^^izlilz£lua^^2 



d-2 I6I74 

1 {1-ny 



±_ln' + 1 ) Fn + (1 - 0^)2^22 - (1 - n')Fu 



+ 



'12 



In 



d-2 4f72 
d-l (1-^^2)2 

'd-2 161^4 

1 l-Jl^ 



Fn 



(65) 



(l-0^)(Fn + F, 



22 i 



n^ + 



d-l 



n' + i]F 



^12 



d-2 402 

1 (1-02)2 



Fo 



(Fii + F22) + 



1 1-0^ 



'12 



d-2 402 ^ ^^ ' ^^^ ' ^_2 '^ d-2 402 
Using the results from sections ^.l] and |5.2| we can express the F's in terms of scalar integrals as, 

Fn = ^Ji(m) + ^i^^J2(A2,m) + ^X3(5^,Ai,A2,m) 

UJ At+CjAi An 

F22 = -Ti(m) + ^ iT2(Ai,m) + ^J3(0,Ai,A2,m) 

F12 = -Ji(m) + -(AiJ2(Ai,m) + A2J2(A2,m))+^-^J3(0,Ai,A2,m) 



(66) 



Fo = m2j3(0,Ai,A2,m) . 

Here we are using a mixed notation, in terms of both lo and O (see (|l3|)), for brevity. Equations 
(63), ( |65[ ) and ( |66[ ) give an explicit analytic expression for 13^'^. Notice also the symmetry of 
(66) under exchange of Ai and A2. The same as with the scalar integral, there are a number of 
particular cases of interest which we briefly comment upon in the remainder of this section. 



5.3.1 The zero recoil case 

In order to study the zero recoil case, it is convenient to write I^^'^ in terms of vectors v'^ 



1/2(^2^ lb V2). Instead of (33) we then have. 



Is'''' = I++vlvK + I^^v'tv'i + I+^v^^v""} + /09''" 



(67) 



v'l = w^. Using the results of sections 5.1 and 5.2 it is not 



In the zero recoil limit, Vi = v[ 

difficult to show that when O = we have v^Xs^^, = 0. Therefore, we can write. 



T,^'^ 



Q=0 



U 



f7=0 



yt^v" + h 



^M^ 



n=o' 



These form factors can be computed as before, resulting in. 



n=o 



n=o 



1 A e \ [^ , , Af-m2 



I^(-i)^'<"'>-^73A; 



2'2(Ai,m) 
2 4m2 



772 



Ai -A. 



■:^2(A2,m)) 



(68) 



(69) 



■J2(Ai,m) 



A2 
^2 



Am' 



Ai -Ao 



■J2(A2,m)) . (70) 
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We notice that we could have arrived at these equations by using a relation analogous to (| 
namely, 



flU 



{li^^'iv^, Ai, m) - T2^'%v'^, A2, m)) 



(71) 



n=o Ai — A2 

showing the consistency of our result. These expressions acquire a particularly simple form for 
some special values of the parameters. For instance, setting m = 0, A^ > 0, j = 1,2 and using 



either (|7lD and (|55|), or (|69[), we get 



h 



-1 



1 



Q=0 
m=0 



/^ 



u 



3 • 1287r2 Ai - A2 
1 1 I 



Af ,..0,^ 






A? 



3^ 
13' 



-H^^-l, 



-M^M. 



3-327r2 Ai- A2 i^ yc y^^^ 

which provides another cross-check of our previous equations. 



^2 



^.oi?; . 



3y 

13^ 
6 



(72) 
,(73) 



5.3.2 The case Ai = = A2 

Another case where the form factors for T^'^^ take a very simple form is when Ai = = A2. In 
this case equations (|6^), ( |65|) and (|66|), together with our previous results for the scalar integrals, 
give. 



h 

hi 

I22 

I12 



m2 i-n"^ 



A,=o 1287r2 n 
2 



log{a^ 



UUK].i 



m^ 



m 



1287r2 

2 
m 

1287r2 



1-^4 (1-^2)3 

+ .^. log(a+) 



2172 



(1-^2)2 
2172 



4f73 



(l-J72)(l-0^ 

4?73 



'^+10,(4' 



log(a+) 






(74) 



We have also computed I^^'^ for A^ = directly from its definition (^) by using the HQET 
method for combining denominators. Full agreement with (^) was found. 

5.3.3 The case m = 

A direct calculation of I^^^ at m = with the HQET method is considerably more involved than 
in the previous case. The results are also much less compact. As an illustration, we will quote the 
result for the form factor Iq when m = and Ai = A2 = A > 0, 



lo 



A2 



2567r2 



<'-"''l^ + ^'°'5<"^ 



hMi 



'1 - Q."^) 
+ 3 + ^ „^ ' log{a+) 



log —— log(a+) + ,^ ( L12 



4n 



2n 



Li2 



2J1 

1 + n 



(75) 



This equation agrees with the general result given by (B^) and (66), evaluated at ttt, = 0, Ai = A2, 
as it should. 
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5.3.4 The chiral logs 

In this section we focus on the case Ai = A2 = Sm, with < Sm/2 < m. We expand in powers 
of 6m/ m, retaining only those terms proportional to log(7n), with coefficients computed to lowest 
order in 6m/ m. In this way, we obtain the chiral logs in I^^'^ . 
From equations (^) and (^3|) , we get, 



Tn I u \ 6m, I u \ 

where the ellipsis denotes terms not containing log(m), or containing higher powers of 6m. On the 
other hand, I3 is needed only to zeroth order in 6m because it enters the form factors with 6m? 
as a coefficient. From (p9|), 

X3(^,0,0,m) = ^i^log(a,)log(^^)+... . (77) 

With these approximations, we obtain the form factors as, 

^11 = ^22 = ^iog(^){-l^-' + ^-'(r^ + ^i°gK))| + --- (78) 

F12 = ^log(5) |-2m2 + 5„^2|^2 + L^log(a,)) K--- (79) 

log p^ m^—— log(a+) + • • • . (80) 




647r2 ^\m'^J n 

These results agree exactly with those of [^, once we take into account the differences in normal- 
ization and conventions. 

6 Final Remarks 

In phenomenological applications, the exact functional dependence of Feynman integrals on masses 
and residual momenta is usually not needed. Often, the first few terms in a series expansion in some 
of the parameters provides the required accuracy. We believe, however, that the exact analytic 
computation presented here does not require more calculational effort than approximate schemes. 
It has the added advantage of being valid over the entire physical region for internal and external 
masses. 

Our result ( |30|) for the scalar three-point integral involves four dilogarithms. This is to be com- 
pared with the analogous vertex integrals in renormalizable theories, which are generally expressed 
in terms of twelve dilogarithms and a collection of logarithms Q, y]. This simplification is afforded, 
of course, by the effective theory formalism, which focuses only on the relevant degrees of freedom. 
Equation (|30| ) is quite compact. Once the values for internal and external masses are given, so 
that the appropriate branches of square roots, logs and dilogs are determined with the aid of the 
"ie" prescription, the expression for T3 given by (Bol) is easily translated into computer code. 

Another possible approach to the computation of three-point integrals is to consider them 
strictly within the context of the effective theory, without introducing auxiliary integrals such as 
T3, eq. (y). In that case, one can parametrize the integrand with the HQET method. The resulting 
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expressions are, however, difficult to handle and, in general, they seem to lead to hypergeometric 
functions of two variables or, more likely, to series of hypergeometric functions. The procedure 
adopted in this paper avoids those difficulties. 
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